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Abstract 

We study QCD with two Dirac fermions in the adjoint representation at finite tem- 
perature by Monte Carlo simulations. In such a theory the deconfinement and chiral 
phase transitions occur at different temperatures. We locate the second order chiral 
transition point at (3 C = 5.624(2) and show that the scaling behaviour of the chiral 
condensate in the vicinity of (3 C is in full agreeement with that of the 3d 0(2) uni- 
versality class, and to a smaller extent comparable to the 3d 0(6) class. From the 
previously determined first order deconfinement transition point (3d = 5.236(3) and 
the two-loop beta function we find the ratio T c /T d « 7.8(2). In the region between 
the two phase transitions we explicitly confirm the quark mass dependence of the 
chiral condensate which is expected due to the existence of Goldstone modes like 
in 3d 0(N) spin models. At the deconfinement transition (ipip) shows a gap, and 
below (3d, it is nearly mass-independent for fixed (3. 
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1 Introduction 



The classification of the chiral transition in finite temperature quantum chromody- 
namics (QCD) for two flavours of quarks is a long-standing problem. This refers 
in particular to analyses of lattice data obtained from simulations with staggered 
fermions. For the continuum theory a second order chiral transition is predicted to 
belong to the universality class of the three-dimensional 0(4) spin model [1, 2, 3]. 
Due to the lattice formulation of the fermionic sector of the theory, chiral sym- 
metry is lost for Wilson fermions and reduced to an 0(2) symmetry for staggered 
fermions. Full symmetry is then recovered only in the continuum limit. Surprisingly, 
corresponding tests with Wilson fermions [4, 5] confirmed the proposed universality 
class, whereas the first tests with staggered fermions and the standard QCD action 
remained unsatisfactory [6]-[10]. Essentially, the comparison of universal critical 
quantities from 0(N) spin models to lattice QCD data failed in the close vicinity 
of the chiral transition. However, at temperatures above the transition some agree- 
ment was found on the pseudocritical line, namely for the form of the line itself and 
for finite size scaling of the chiral condensate on the line [11]. The apparent failure 
of the tests was attributed to a very small critical region, i. e. use of too large m q a- 
values and lattice artifacts, i. e. calculations with too small N T and/or unimproved 
actions. 

The issue is still under consideration and still unsettled: Kogut and Sinclair, who 
use a QCD action including an irrelevant 4-fermion term, conclude from simulations 
at m q = on a 24 3 x 8 lattice [12] that the critical region must be very small, 
but the transition class could be compatible to the 3d 0(2)/ 0(4) universality class. 
Chandrasekharan and Strouthos [13] performed simulations at strong coupling and 
compared their data to predictions from chiral perturbation theory for 0(2). They 
find full agreement for masses which are much smaller than those normally used in 
QCD calculations. Contrary to that, D' Elia et al. [14] exclude 0(2)/0(4) scaling 
in their finite-size-scaling analysis of specific heat and susceptibility data for the 
standard QCD action. Instead they find consistency of their data with a first order 
transition. 

The difficulties and discrepancies encountered in the investigation of the nature 
of the chiral transition could be related to another puzzle of finite temperature QCD: 
the coincidence of the chiral and the deconfinement transitions. Though there is no a 
priori reason why the two transitions should occur at the same temperature in QCD, 
they seem to be strongly bound to each other. The interplay of the deconfinement 
and chiral transitions could then obscure the critical chiral behaviour in the vicinity 
of the common critical point. It is possible to study the properties of the chiral 
transition without direct interference from the deconfinement transition. In a special 
QCD-related model, SU(N C ) gauge theory with two Dirac fermions in the adjoint 
representation (aQCD), the deconfinement transition occurs at a temperature 
which is definitely smaller than T c , the chiral transition temperature [15, 16]. 
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In general, aQCD with Nf massless fermions exhibits an SU (2Nf) chiral symme- 
try, which is spontaneously broken at low temperatures to SO(2Nf) [17, 18]. Since, 
in contrast to QCD, the fermionic part of the aQCD action does not break the global 
Z(N C ) centre symmetry of the SU(N C ) gauge part, the model holds in addition the 
same global Z(N C ) symmetry as pure SU(N C ) gauge theory. Correspondingly, one 
expects two finite temperature transitions and there are two independent order pa- 
rameters, the Polyakov loop and the chiral condensate, which indicate the break- 
ing of the respective symmetries. The universality class of the chiral transition of 
aQCD is normally different from that of QCD. For two flavours of adjoint Dirac 
fermions the chiral symmetry group is SU(4). That group is isomorphic to 0(6), 
which breaks to 0(5), implying the existence of 5 Goldstone bosons. The breaking 
SU(A) — > SO (A) leads however to 9 Goldstone bosons. The corresponding transi- 
tion belongs therefore to a different universality class. In a recent paper [19] this 
issue has been investigated using renormalization-group arguments and first results 
for the critical exponents of this class have been obtained: rj ps 0.2 , v ps 1.1. The 
points we have just raised pertain to continuum theory. On the lattice, the stag- 
gered formulation leaves, like in QCD, only a remaining 0(2) chiral symmetry. In 
the continuum limit one should however recover the full SU(A) symmetry and the 
critical behaviour following from its breaking to 5*0(4). 

The main objective of this paper is to investigate the chiral transition of aQCD 
and if possible to establish its universality class. Since we are simulating aQCD 
on the lattice with two flavours of staggered Dirac fermions 1 we expect to find the 
three-dimensional 0(2) class. The thermodynamics of aQCD with N c = 3 has been 
studied to a great extent already by Karsch and Liitgemeier [16]. In particular, 
a first order deconfinement transition was found at Pa = 5.236(3) and a second 
order chiral transition at /3 C = 5.79(5) was inferred from the mass dependence of the 
chiral condensate and its susceptibility. However, no comparison to the universal 
scaling functions of 0(N) spin models was attempted in the vicinity of f3 c . In order 
to achieve that we have made additional simulations at more /3-values around the 
supposed chiral transition point and at additional and smaller m^a-values. A second 
aim of our paper is the study of the Goldstone effects which are connected to the 
breaking of chiral symmetry. If aQCD behaves as an effective three-dimensional 
0(N) spin model close to the chiral transition point then one would expect the 
chiral susceptibility Xm to diverge below (3 C as (m^a) -1 / 2 for m q a — > 0. Such a 
mass dependence has in fact been observed in the range (3d < (3 < f3 c in Ref. [16]. 
For smaller T- values, approaching T = 0, the theory changes to an effective four- 
dimensional theory and the expected divergence of Xm becomes logarithmic in m q a. 
The change in the behaviours of (ipip) and Xm with decreasing T, also below (3d, is 
therefore of interest. 

1 We have adopted the QCD convention of counting the number of fermion flavours. Because 
the fermions are in the adjoint representation one Dirac fermion corresponds to two Majorana 
fermions. 
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The paper is organized as follows. In the next section we define the action and 
the observables of aQCD, which we want to consider. Then we discuss numerical 
details of our simulations and give an overview of the available data. In Section 3 
we recapitulate those results from the O(N) spin models which are most relevant for 
the subsequent analysis of the aQCD data in Section 4. We close with a summary 
and the conclusions. 



2 QCD with adjoint fermions on the lattice 

In the following we shall use previous results obtained by Karsch and Lutgemeier 
and extend their simulations for our purposes. The version of aQCD which we use 
is therefore the same as the one already described in detail in Ref. [16] and it will 
suffice to address only the main aspects of the corresponding lattice theory. The 
action is similar to QCD with staggered fermions in the fundamental representation. 
The gluon part S G of the action is the standard Wilson one-plaquette gauge action 
with SU(3) link variables and coupling (3 = Q/g 2 

S G = pJ2{ 1 -l RetTlJ p ) } ■ C 1 ) 
p 6 

In contrast to QCD the fermions are in the eight-dimensional adjoint representation 
of the SU(3) group and have eight colour degrees of freedom. The fermionic part 
Sp of the action requires then also an eight-dimensional representation of the 
link matrices which may be expressed in terms of the three-dimensional matrices 

U$ = ltr[X a U^X b U^} , (2) 

where the A a are the eight Gell-Mann matrices and tr is as in (1) the 3-trace. The 
total aQCD action is 

S = S G + S F = S G + J2^M(U^) x ^ y . (3) 

Here, M(U^) is the standard staggered fermion matrix where the link variables 
have been replaced by the respective matrices. We note some important 
properties of . From Eq. (2) it follows that is real and its 8-trace (tr 8 ) is 
related to the 3-trace of via 

tr 8 U^ = \trU^\ 2 -l . (4) 

Further, is invariant under Z(3) transformations of the U^> -matrices because 
the centre elements appear in Eq. (2) as a complex conjugate pair. As a consequence, 
Sf does not break the Z(3) centre symmetry and the Polyakov loop 
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X X4=l 



(5) 
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is an order parameter for the deconflnement transition like in pure SU(3) gauge 
theory in the limit N a — > oo. The order parameter for the chiral transition, the 
chiral condensate, is defined by 

<^> = -^ = ^^<TrM- 1 >. (6) 
XV ^ 7 V dm q a N*N T 2 V 1 y ! 

where the trace Tr has to be taken over the eight colours and all spacepoints. The 
derivative of the chiral condensate with respect to m q a is the chiral susceptibility 

diW) _ T d 2 \nZ 
^ m dm q a V d{m q a) 2 

It is the sum of the connected term Xcon and the disconnected term Xdis- The latter 
measures the fluctuations of the chiral condensate. The two parts are given by 

Xais = j^(^)\((^M-r)-(TrM- i r) , (8) 

Xcon = -^J2( M *,0 M °1) ■ ( 9 ) 

x 

In our subsequent analysis we shall utilise the data for the Polykov loop, the chiral 
condensate and the disconnected part of the susceptibility. 

Since the matrices are real this is also the case for the total fermion matrix 
M(U^). The determinant of M(U^) may then be represented as a Gaussian 
integral over real bosonic fields with the advantage that the exact Hybrid Monte 
Carlo algorithm can be used for the simulation of aQCD with two adjoint fermion 
flavours. The corresponding effective pseudo-fermion action is given by 

S eff = S G (U i3) ) + ¥ (M(U i8) YM(U {8) )) $ , (10) 

with real eight-component <3>-fields, which reside only on even lattice sites to com- 
pensate for the fermion doubling introduced by the use of M l M instead of M. The 
molecular dynamics equations for aQCD are more complicated as the ones for QCD, 
because M is now quadratic in U^. Details of the update algorithm can be found 
in appendix A of Ref . [20] . 

Our simulations were performed on N% x 4 lattices with N a = 8, 12 and 16 and a 
fixed length r = n5r = 0.25 of the trajectories. The number of molecular dynamics 
steps n was chosen such that the final Metropolis acceptance rate was in the range 
of 60-80%. On the 8 3 x 4 lattice n varied between 20 and 40 for 0.10 > m q a > 0.02. 
For lower masses and/or increasing spatial extension N a a of the lattice the number of 
steps had to be increased. The optimal n-value was however relatively independent 
of /3, apart from the region around and below the deconflnement transition, where 
larger n were necessary. The chiral condensate and the Polyakov loop were measured 
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Figure 1: The chiral condensate (ipip) (a) and the Polyakov loop L3 (b) versus (3 
and m q a on the 8 3 x 4 lattice (plot (b) has been rotated by 90°). The filled squares 
denote the new data, the empty ones are from Ref. [16]. The dotted lines indicate 
the height. The data points are connected by straight lines to guide the eye. 



5 



after the completion of each trajectory. Both the condensate and the disconnected 
part of the susceptibility were calculated from a noisy estimator [21] with 25 random 
source vectors. As already mentioned in the introduction, our first intention was 
to enlarge the data set of Ref. [16]. To this end we have performed additional 
simulations on the 8 3 x 4 lattice with 900-2000 trajectories for measurements at 
m q a = 0.005,0.03 and 0.06, and for the new /9-values 5.1, 5.45, 5.55, 5.6 and 5.65. 
The already existing data at m q a = 0.01, 0.02, 0.04, 0.05, 0.08 and 0.10 and at the 
/^-values 5.2, 5.3, 5.35, 5.4, 5.5, 5.7, 5.75, 5.85, 5.9, 6.2 and 6.5 have been completed 
to obtain a mesh covering the region from above the chiral phase transition to a 
point just below the deconfinement transition. In order to see possible finite size 
dependencies we have carried out further simulations at selected points on larger 
lattices with N a = 12 and 16. The respective data will be show in the figures of 
Section 4. The chiral condensate data from the different lattices coincide inside the 
error bars at essentially all points. In the chiral susceptibility data there are however 
some differences which may be due to insufficient statistics. 

An overview of the data from the 8 3 x 4 lattice is given in Fig. 1 where we show 
three-dimensional plots of the chiral condensate (ipip) and the Polyakov loop L3. In 
the plots we have indicated on the line at m q a = the previously [16] determined 
transition points (3d = 5.236(3) and (3 C = 5.79(5) by empty squares. The filled 
square on the same line denotes our new estimate for (3 C , which we will explain in 
Section 4. Obviously, the chiral phase transition point cannot be determined by 
a simple inspection of the order parameter (ipip), even though we have many new 
data points in the critical region. For the first order deconfinement transition the 
situation is much clearer. The Polyakov loop L 3 shows jumps between (3 = 5.2 and 
5.3 , depending on the mass value. Our new results at (3 — 5.1 and 5.2 are definitely 
below the deconfinement transition. Insofar they support the determination of the 
exact position by Karsch and Liitgemeier who produced more data in the transition 
region at the masses m q a = 0.02, 0.04, 0.08 and 0.10 . The straight line connections 
in Fig. 1 in this region at the other mass values should therefore be taken with care. 
Two further comments are appropriate here. Already in Ref. [20] a remarkable fact 
was noted: the first order deconfinement transition point (3d = 5.236(3) of aQCD 
coincides with the second order chiral transition point (3 C = 5.233(7) of QCD with 
two fundamental quarks and the standard staggered action [9]. In addition, we see 
in Fig. 1 that {^ip) exhibits also a discontinuity at (3d , whereas the Polyakov loop 
shows no sign of the chiral transition. 



3 The three-dimensional 0(N) spin models 

In this section we want to summarize the basic facts about the three-dimensional 
O(N) spin models and their universality classes. A more general survey on O(N)- 
symmetric systems can be found in Ref. [22]. The numerical results to which we will 
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refer in the following have been obtained from simulations of the standard O(N)- 
invariant non-linear sigma models with the partition function 



cxp 



j (fx ■ 4> y + h ■ <p x 



(11) 



Here, x and y are nearest-neighbour sites on a three-dimensional hypercubic lat- 
tice and <f) x is an iV-component unit vector at site x. The coupling J acts as inverse 
temperature J = 1/T and H is the magnetic field. The spin vector <p x may be decom- 
posed into a longitudinal (parallel to H) and a transverse component (perpendicular 

to H): <f) x = (f>xH + <j) x . The order parameter of the system, the magnetization M 
and its susceptibility xl are then defined by 

M = U £ *i) = <<*"> • XL = |f = V - W) , (12) 

that is, M is the expectation value of the lattice average of the longitudinal spin 
component and xl its derivative with respect to H — \H\ (implying H > 0) . 

The spin system exhibits a second order phase transition at a critical temperature 
T c = 1/J C . For H — > and T < T c the 0(iV)-symmetry of the spin system is 
spontaneously broken and the magnetization attains a finite value M(T, 0) > 0, 
whereas M(T > T c , 0) = 0. The behaviour of M(T, H) in the region around the 
critical point is described by the equation of state 

M = h 1/s f G {z) , with z = t/h 1/f3S , (13) 
where /g(z) is a universal scaling function, which is normalized such that 

/ G (0) = 1 , and f G (z) zZ ^J-zf . (14) 
This is achieved by the introduction of the reduced temperature t and field h with 

t = (T- T c )/T and h = H/H , (15) 

such that 

M(t = 0, h) = h 1/s , M(t < 0, H -> 0) = (-if . (16) 

The critical exponents 5 and /5 determine all other critical exponents via the hyper- 
scaling relations (here d — 3) 

7 = /3(<J-1), di/ = + = 2 -a, 2-77 = 7 A' ■ (17) 

The scaling function /g(^) and the critical exponents are the same for all members 
of the corresponding universality class. They can therefore be used in tests on the 
class of a second order phase transition. In a series of papers these functions have 
been calculated for the three-dimensional O(N) models with N = 2 [23], N = A 
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Figure 2: The universal scaling function fo = Mh l l 5 versus the scaling variable 
z = th~ 1/f3S for N = 2 [23], N = A [26] and iV = 6 [27]. Parametrizations of the 
curves can be found in Refs. [23, 26, 27]. 

[24, 25, 26] and N = 6 [27]. They are shown in Fig. 2. We observe similar shapes 
of the different /c(^)-curves with an increasing slope for larger N. In the plot the 
critical line T = T c is mapped to the point z — 0, the symmetric (high temperature) 
phase T > T c , H — > corresponds to the limit z — > oo and the coexistence line 
T < T c , H — ► to 2 — > — oo. Since the exponents 5 and /? are slightly different for 
each N the scales of both axes of the plot vary with N. In Table 1 we have listed 
the critical exponents which have been used in Fig. 2. We see that the exponent /3 



N 


P 


5 


V 


a 


1/5 




Ref. 


2 


0.349 


4.7798 


0.6724 


-0.017 


0.2092 


0.5995 


[23] 


4 


0.380 


4.824 


0.7377 


-0.213 


0.2073 


0.5455 


[26] 


6 


0.425 


4.780 


0.818 


-0.454 


0.2092 


0.4922 


[27] 



Table 1: The critical exponents of the 3d O(N) spin models. 
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is increasing with N, while 5 is nearly constant. The main influence of N in a 
scaling plot is therefore coming from the iV-dependence of the product A = (35 
which appears in the scaling variable z. 

Since the susceptibility xl is the derivative of M with respect to H, it can as 
well be expressed in terms of a universal scaling function 

Q M h 1 ' 8 - 1 . . 

XL = ~dW = lir^ ■ (18) 

with 

/*(*) = ] (fa{z) ~ • (19) 

The scaling function f x (z) has a peak at a positive value z = z pc . Correspondingly, 
Xl has a maximum for each fixed value h at a temperature t pc (h). This function of 
h defines the pseudocritical line in the (t, /i) -plane for t > 0. Obviously, it is given 
by 

t pc /T 1/A = v or T pc -T c = cH 1 ^ , (20) 

where c = ZpcT^H^ 1 ^ is a non-universal constant. Eq. (20) is often used to deter- 
mine T r . 



Due to the existence of massless Goldstone modes in O(N) models with N > 1 
and dimension 2 < d < 4 singularities are expected on the whole coexistence line 
T < T c , H = [28, 29, 30] in addition to the known critical behaviour at T c . In 
particular, the susceptibility xl is diverging on the coexistence line. For all T < T c 
and H — > the divergence is 

N f #- 1/2 for d = 3 
X L(T,iJ)~ , (21) 

! In if for d — 4 

which implies the following //-dependence at fixed T for the magnetization 

/ x / x f # 1/2 for rf =3 

M(T,/J)=M(T,0)+ Cl (T). . (22) 

I H{mH — 1) tor d =4 

For the three-dimensional O(N) models with = 2,4,6 the i/ 1//2 -dependence of 
the magnetization for T < T c has been verified numerically in Refs. [23, 25, 27]. In 
all cases one finds that the slope Ci(T) is small for very low temperatures, that is 
M is essentially constant. Approaching T c , the slope increases and the region close 
to H = where M is a straight line as a function of if 1 / 2 shrinks. This is of course 
what one expects, because at T c the magnetization changes its if -dependence from 
H 1 / 2 to H l / & with the effect that c\ — > oo. The scaling functions f G and f x describe 
this behaviour correctly in the critical region. Yet, Eqs. (21) and (22) are also valid 
for all temperatures below the critical region. 
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4 Analysis of the aQCD data 



In aQCD the temperature is given by T = 1/ ( N T a) and the analogon of the magnetic 
field is the quark mass m q a. Since the lattice spacing a is a monotone function of the 
coupling (3, we can use f3 instead of the temperature and temperature differences are 
proportional to coupling differences. Accordingly, we define instead of the reduced 
temperature t a reduced coupling (3 r and a reduced field m, where 

(3 r — {(3 — /3 c )//3q , m = (m q a) / (m q a) and z = fl r m~ x ^ . (23) 

The constants (3q and (m q a)o have to be fixed from the normalization of the scaling 
function for (iptp). It is clear, that before any scaling tests can be carried out, 
the critical point (3 C of the chiral transition has to be determined as accurately as 
possible. 

4.1 Location of the chiral transition point 

We start our search by following first the usual strategy to extrapolate the pseudo- 
critical couplings (3 pc (m q a) to m q a = with an ansatz equivalent to Eq. (20) 

(3 pc = (3 c + c(m q a) 1 ' A . (24) 

The (3 pc are in principle the positions of the maxima of Xm- I n Ref. [7] the two parts 
Xdis and Xcon of Xm have been studied seperately. It was found that the contribution 
of the connected part to Xm is a slowly varying function of (3 in the critical region and 
that the ratio Xcon/ Xdis decreases with decreasing quark mass. It suffices therefore to 
take the peak positions of Xdis for the extrapolation to the critical point. The data 
for the susceptibility Xdis from the 8 3 x 4 lattice are shown in Fig. 3 for the three 
smallest masses 0.005, 0.01 and 0.02. They have been interpolated by reweighting 
where we had enough statistics. We observe a pronounced but broad maximum for 
the smallest mass 0.005 around (3 rs 5.73, for m q a = 0.01 a less peaked and smaller 
maximum around (3 ~ 5.74 and for m q a = 0.02 only a long and low plateau between 
(3 = 5.65 and (3 = 5.85 . In Table 2 we have listed our estimates for the pseudocritical 
couplings. For m q a = 0.01 and 0.02 they differ somewhat from the values obtained 
in Ref. [16] because we have more data points, the value at m q a = 0.005 is new. 
Certainly, the errors on (3 pc (m q a) are too large to lead to a precise value of the 



m q a 


0.005 


0.010 


0.020 


(3pc 


5.73(8) 


5.74(6) 


5.75(10) 



Table 2: The pseudocritical couplings (3 pc from the 8 3 x 4 lattice. 
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Figure 3: The susceptibility Xdis as a function of (5 for the three smallest masses 
m q a = 0.005, 0.01 and 0.02 on the 8 3 x 4 lattice. The solid lines are the results from 
reweighting, the dashed lines show the Jackknife error corridor. 



critical point j3 c . On the other hand it is clear that (3 C must be definitely smaller 
than the (3 pc {m q a) because the susceptibility peaks only in the high temperature 
phase (3 T > 0. If one nevertheless tries an extrapolation with Eq. (20) using the 
numbers from Table 2 and 0(2) and 0(6) values for A one obtains (3 C = 5.72(16) 
and 5.71(19), respectively. 

A second possibility to locate the critical point is based on the expected mass 
dependence of the order parameter at (3 — (3 C 

(^) = d c {m q a) l ' & . (25) 

The method consists in a comparison of the mass dependence of the data at each 
/?-value with this functional form. In practice one makes straight line fits at fixed j3 
of the logarithm of the order parameter as a function of ln(m 9 a) 

In (V>V) = In d c -\- — ln(m„a) . (26) 
o 

In the vicinity of f3 c one should find the lowest x 2 '/d.o.f. and/or the largest fit range 
in m q a. In addition, one obtains a first result for 1/5. A method of this type was 
successfully applied to the corresponding finite size scaling problem in SU (2) gauge 
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p 


In d c 


1/5 


m q a-i ange 


N p 


Xminl 1 d.O.f '. 


Xmaxl d.O.f. 


5.80 


1.510(37) 


0.4246(97) 


0.005 - 0.03 


4 


7.51 


32.63 


5.75 


1.336(25) 


0.3549(65) 


0.005 - 0.03 


4 


4.44 


28.00 


5.70 


1.124(23) 


0.2786(60) 


0.005 - 0.03 


4 


1.05 


10.43 


5.65 


1.026(09) 


0.2322(26) 


0.005 - 0.06 


7 


1.64 


3.00 


5.60 


0.949(04) 


0.1921(17) 


0.005 - 0.10 


9 


5.25 


9.11 


5.55 


0.822(19) 


0.1423(50) 


0.005 - 0.04 


5 


4.18 


7.51 


5.50 


0.755(19) 


0.1083(46) 


0.005 - 0.03 


4 


5.62 


8.79 



Table 3: The parameters of the fits according to Eq. (26). The number N p of fitted 
data points and the fit range correspond to the minimal x 2 /d.o.f. with N p > 4. The 
last column is the maximal x 2 1 'd.o.f '. , if N p is varied between 4 and 9. 



theory [31]. In Table 3 we present the results for all fits in the interval 5.5 < (3 < 5.8. 
At each /5-value we have nine data points in the range 0.005 < m q a < 0.10. Because 
we have two fit parameters we used at least the data points from the four smallest 
masses, if possible more, up to all the nine points. The best x 2 1 d.o.f . is obtained 
at f3 — 5.7 and 5.65, but for f3 = 5.7 only in the smallest mass interval. At (5 — 5.65 
there is only a small variation in x 2 1 d.o.f. when N p is increased. For f3 = 5.60 we 
find a larger x 2 /d.o.f., which is due to the lowest mass point, but the best fit includes 
all nine data points. The data points at m q a = 0.005, 0.01 and 0.02 from the 8 3 x 4 
lattice agree here within the errorbars with those from the 16 3 x 4 lattice, so that 
a finite size effect is improbable. In conclusion, the range 5.6 < j3 < 5.7 contains 
presumably the critical point f3 c . If we look at the corresponding 1/5- values from 
the fits we find 0.1921 < 1/5 < 0.2786. That interval includes the value 0.2092, 
which is expected from the 0(2) and 0(6) spin models, suggesting 5.6 < f3 c < 5.65. 

We can further improve the search for the critical point by assuming a scaling 
ansatz for the order parameter 

<^> = m l,5 f{z) (27) 

as in spin models. Expanding f(z) for small \z\ one obtains 

= m 1 / 5 (/(0) + ^'(0) + ...) (28) 

= ( mq a) 1 l & (d c + d\{f3-f3 c ){m q a)- l l A + ..) . (29) 

Here we have used the normalization condition /(0) = 1 and Eqs. (23) and (25). 
They imply d c = (m q a) and a similar relation for the amplitude d\ . The ad- 
vantage of Eq. (29) is, that it is valid also in the vicinity of f3 c and at medium size 
values of m q a. The mass should not be too small to keep \z\ small and not too 
large to be inside the critical region where (27) is applicable. We have made fits 
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p 


d c 




m q a-i ange 


X 2 /d.o.f. 


5.90 


2.736 (8) 


-0.1101 (16) 


0.03 - 0.10 


2.19 


5.85 


2.704 (8) 


-0.0868 (14) 


0.02 - 0.10 


1.53 


5.80 


2.706 (6) 


-0.0713 (12) 


0.02 - 0.10 


1.16 


5.75 


2.689 (6) 


-0.0508 (11) 


0.02 - 0.10 


1.30 


5.70 


2.662 (6) 


-0.0283 (10) 


0.01 - 0.10 


1.56 


5.65 


2.654 (5) 


-0.0089 (08) 


0.005 - 0.10 


1.68 


5.60 


2.659 (7) 


0.0080 (13) 


0.01 - 0.10 


2.97 


5.55 


2.624 (9) 


0.0342 (20) 


0.03 - 0.10 


2.30 


5.50 


2.606(13) 


0.0583 (24) 


0.03 - 0.10 


1.99 



Table 4: The parameters for the fit (29) with 0(2) exponents. 



p 


d c 


dl(P-Pc) 


m^a-range 


X 2 /d.o.f. 


5.90 


2.862 (8) 


-0.1722 (16) 


0.02 - 0.10 


2.43 


5.85 


2.811 (9) 


-0.1449 (21) 


0.01 - 0.10 


2.89 


5.80 


2.793 (6) 


-0.1192 (17) 


0.01 - 0.10 


0.75 


5.75 


2.748 (6) 


-0.0841 (15) 


0.01 - 0.10 


1.20 


5.70 


2.693 (6) 


-0.0463 (14) 


0.005 - 0.10 


2.39 


5.65 


2.669 (6) 


-0.0159 (15) 


0.005 - 0.10 


1.73 


5.60 


2.647 (8) 


0.0141 (22) 


0.01 - 0.10 


2.70 


5.55 


2.585(11) 


0.0563 (34) 


0.03 - 0.10 


2.07 


5.50 


2.538(15) 


0.0961 (40) 


0.03 - 0.10 


1.48 



Table 5: The parameters for the fit (29) with 0(6) exponents. 



with Eq. (29) and 0(2) and 0(6) exponents alternatively. The results are given in 
Tables 4 and 5. The shown fit range in m q a extends always to m q a = 0.10. The 
lower limit of the range was chosen such that the next smaller mass would increase 
X 2 /d.o.f. considerably. We observe that the ranges extend to the smallest masses 
for 5.6 < P < 5.7. The result for the amplitude d c should of course be independent 
of the used (5- value as long as (5 is in the critical region. For 0(2) exponents (Table 
4) the amplitude result is generally less varying than for 0(6) exponents (Table 5), 
in the range 5.6 < j3 < 5.7 the result for 0(2) is actually constant within error 
bars. If we consider this /5-range as the effective critical region we obtain the value 
d c = 2.66(1) for 0(2) and d c = 2.67(3) for 0(6) exponents. The assessment of the 
range 5.6 < (3 < 5.7 as the critical region is further supported by the parameter 
d\{(5 — p c ), because it has a zero in this region. In Fig. 4 we have plotted the results 
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Figure 4: The parameter d\{(3 — (3 C ) as a function of (3 for 0(2) (circles) and 0(6) 
exponents (triangles). The straight lines are fits to the values at (3 — 5.6,5.65 and 
5.7 . The zeros of the fits are given by the numbers in the plot. 



for this parameter and straight line fits to the three values in the critical region. We 
obtain coinciding values for (3 C : 5.6240(23) for the 0(2) and 5.6235(2) for the 0(6) 
case. The errors on these results should however not be overestimated, because 
they are only those from the straight line fits. Nonetheless we have arrived at a 
remarkably clear estimate for f3 c . From f3 c and (3d we can derive the ratio T c /T d of 
the critical temperatures approximately with the 2-loop /3-function 



aA L *R W = { ^ ^[~). (30) 



The coefficients bo and b\ for aQCD are given by [32] 

3 b 90 

167T 2 ' 1 (167T 

and the temperature ratio becomes 



b = 77~5. > fe i = ~ n^2\2 > ( 31 ) 



exp —{(3 C - (3 d ) ) = 7.8(2) . (32) 



T d a((3 c ) \(3 d J " V 9 

Though we have a lower estimate for (3 C than Ref. [16] we nevertheless arrive at 
about the same ratio. That is due to the different value for b±, which was wrongly 
taken as +115/(384tt 4 ) in [16]. 
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4.2 The Goldstone effect 



We have just successfully applied the scaling ansatz (27) to the chiral condensate 
data. One might therefore directly try to compare the data with the universal scaling 
functions fc{ z ) of the spin models. Before we can do so, however, we have to fix the 
normalization f3 of the reduced temperature (3 T (the normalization of m is already 
known from d c ). The amplitude d\ does not fix (3q because the expression zdf /dz in 
Eq. (28) is independent of any normalization of z or j3 r . For that purpose the critical 
amplitude of the chiral condensate on the coexistence line has to be determined. 
In the case of the three-dimensional O(N) spin models the magnetization data for 
T <T C and H > have been extrapolated to H = using Eq. (22). We can proceed 
in the same manner to obtain (ipip) -values at m q a = 0, supposed the corresponding 
Goldstone ansatz 

(W)((3,m q a) = (W)(J3,0) + c^im.a) 1 ' 2 + c 2 ((3)(m q a) + ... (33) 

for /3 < (3 c is meaningful at all. In Fig. 5 we have plotted the chiral condensate data 
for all /^-values between 5.3 and 5.9 versus (m q a) 1 ' 2 . Here, also the available data 
from the 12 3 x 4 and 16 3 x 4 lattices are shown. Since there are no perceptible finite 




0.1 0.2 0.3 



Figure 5: The chiral condensate (ipip) as a function of (m^a) 1 / 2 for all (5- values 
between 5.3 (highest values) and 5.9 (lowest values) from iV 3 x 4 lattices with N a = 8 
(circles), 12 (triangles) and 16 (diamonds). The lines are fits with the ansatz (33), 
the filled circles denote the extrapolations to m q a = 0. 
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p 


(W(/5,0) 


Cl 


c 2 


m g a-range 


Xminl 1 d.O.f '. 


5.30 


1.569 (13) 


1.23 (12) 


- 


0.005 - 0.03 


0.24 


5.35 


1.447 (07) 


1.49 (04) 


- 


0.005 - 0.08 


1.33 


5.40 


1.295 (08) 


1.85 (04) 


- 


0.005 - 0.08 


0.63 


5.45 


1.141 (11) 


2.64 (11) 


-1.72 (25) 


0.005 - 0.10 


3.03 


5.50 


1.020 (21) 


2.81 (26) 


-1.37 (74) 


0.005 - 0.06 


1.19 


5.55 


0.854 (22) 


3.49 (22) 


-2.48 (54) 


0.005 - 0.08 


2.18 


5.60 


0.697 (24) 


3.98 (29) 


-2.66 (84) 


0.005 - 0.06 


1.67 


5.65 


0.468 (17) 


5.48 (20) 


-5.97 (54) 


0.005 - 0.08 


1.33 


5.70 


0.285 (19) 


6.44 (22) 


-7.83 (61) 


0.005 - 0.06 


1.28 


5.75 


0.075 (22) 


7.74 (30) 


-10.77 (96) 


0.005 - 0.05 


0.61 


5.80 


-0.025 (19) 


7.51 (19) 


-8.71 (45) 


0.005 - 0.08 


1.26 


5.85 


-0.190 (12) 


8.39 (14) 


-10.34 (37) 


0.005 - 0.08 


1.14 


5.90 


-0.267 (10) 


8.27 (12) 


-9.46 (31) 


0.005 - 0.08 


1.97 



Table 6: The parameters of the fits according to the Goldstone ansatz (33). The 
upper fit boundary was fixed by the minimal x 2 jd.o.f. given in the last column. 



size effects, the correlation length must still be small even at the lowest masses. 
Because we have no data for masses below m q a = 0.005 we cannot see directly 
from Fig. 5 at which (3 the Goldstone behaviour, Eq. (33), changes to Eq. (25), the 
behaviour at the critical point. Nevertheless we have made fits to the data at each 
f3- value according to Eq. (33), including the first three terms. In Table 6 we present 
the results of the fits, the corresponding curves are shown in Fig. 5. For (3 = 5.3, 
5.35 and 5.4 the parameter c 2 is zero within error bars. At f3 — 5.3 the influence 
of the nearby deconfinement transition can explicitly be seen at the higher masses 
m q a > 0.04 , at (3 — 5.35 for m q a = 0.1. The fitcurves for (3 < (3 C are shown as solid 
lines. The short horizontal line in Fig. 5 separates the extrapolated points below 
and above the critical point. Obviously, the data meet the Goldstone expectations 
for a three-dimensional theory quite well. 

The fitcurves for (3 > (3 C are shown as dashed lines and since they reproduce 
the data for m q a > 0.005 as well with a small x 2 /d-o.f., the fits as such provide no 
clear sign for the location of f3 c . However, they provide an upper limit for (3 C insofar 
as (ipip) at m q a = should be positive. From Table 6 we see that this is not the 
case for (3 > 5.8 . At /3 ~ 5.8 the extrapolation is zero and because of that (3 C was 
located near to 5.8 in Ref. [16]. A hint to the genuine critical point is provided by 
the /9-dependence of the parameter c 2 . With increasing coupling c 2 increases slowly, 
at (3 — 5.65 however it jumps suddenly to a higher value indicating the change of 
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the behaviour from Eq. (33) to Eq. (25). We have also tried fits for 5.3 < (3 < (3 C 
with a Goldstone ansatz corresponding to a four- dimensional theory 

(^){(3,m q a) = (^){(3,0) + d 1 {(3)[{m q a)ln{m q a)-{m q a)} + ... . (34) 

The data are not well reproduced with such an ansatz and the respective y 2 /d.o.f. - 
values are about 10-50 times larger than those in Table 6. That is we confirm that 
aQCD at finite temperatures in the range < T < T c behaves as an effective 
three-dimensional spin model. 

In Fig. 6 we have plotted the disconnected part of the chiral susceptibility Xdis 
for (3 = 5.6 and 5.65, that is just below and above the critical point (3 C . A clear 
finite size effect is not visible, possibly because our statistics are too low. However, 
we notice a change in the mass dependence. Below the critical point Xdis seems to 
follow the behaviour expected for Xm from the Goldstone ansatz (33) 

Xm((3, m q a) =c 2 + y (m g a)" 1/2 , (35) 

whereas above (3 C the increase in the susceptibility as a function of the mass seems 
to be stronger. A similar change in the mass dependence of Xdis from (3 = 5.5 to 
5.7 can be seen in Fig. 4 of Ref. [16]. For comparison we show in Fig. 6 also the 
predictions for Xm following from Eq. (35) and the parameters from Table 6. At 
(3 — 5.6 the agreement with the data suggests that here Xcon is negligible. 

We have also investigated the Goldstone effect at (3 = 5.1 and 5.2 . Both values 
are below the deconfinement transition point (3d = 5.236(3) . In Fig. 7 we show the 




5 10 15 5 10 15 



Figure 6: The chiral susceptibility Xdis versus (m q a) 1//2 for (3 = 5.6 (a) and (3 = 5.65 
(b). The straight lines correspond to Eq. (35) with the parameters from Table 6. 
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Figure 7: The chiral condensate (ipip) for [3 — 5.1 (upper points) and (3 = 5.2 (lower 
points) versus {m q a) l l 2 (a) and versus m q a(ln(m q a) — 1) (b). The straight lines are 
fits to ansatz (33) in (a) and to ansatz (34) in (b). 



data for the chiral condensate. Due to the nearby deconfinement transition the data 
for (3 = 5.2 are still too noisy to enable an analysis. At (3 — 5.1, on the other hand, 
the data are nearly independent of the mass and fits to ansatz (33), Fig. 7(a), and 
ansatz (34), Fig. 7(b), are both possible, without favouring one of them. 

As mentioned already, the verification of the Goldstone ansatz (33) leads to an 
estimate of the chiral condensate on the coexistence line. In the vicinity of the 
critical point we expect the following critical behaviour 

< /3c, m q a = 0) = B((3 C - (3)^ = {-f3 r f™ . (36) 

Here, (3 m is the magnetic critical exponent (called (3 in the case of the O(N) spin 
models). From the second part of Eq. (36) it is clear that the critical amplitude B 
and the normalization (3q are related by (3q = B" 1 '^" 1 . In order to determine B we 
have to evaluate the estimates of the chiral condensate at m q a = from Table 6. 
These data are shown in Fig. 8. Of course, at larger values of f3 c — (3 corrections 
to the leading term are expected to contribute. We have therefore first tried fits to 
Eq. (36) without the last three points, using [3 m from the 0(2) and 0(6) models, 
or as a free parameter. In all three cases the x 2 /d.o.f. is of the order 16-20 when 
the point at (3 = 5.6 is included in the fit. Most probably, due to the vicinity of the 
critical point, the extrapolation of the chiral condensate to m q a — at (3 = 5.6 is 
not as reliable as at the higher (3- values. Discarding that point, we obtain excellent 
and essentially coinciding fits, if [3 m is treated as a free parameter and for (3 m {0{2)). 
The fit is not as good, if (3 m from the 0(6) model is used. We have also made fits 
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B 






9 


(P c - /?)-range 


X 2 /d.o.f. 


2.10(02) 


0(2) 






0.074 - 0.174 


0.075 


2.43(02) 


0(6) 






0.074 - 0.174 


4.100 


2.06(12) 


0.338(30) 






0.074 - 0.174 


0.012 


2.04(06) 


0(2) 


1.39(86) 


1.98(68) 


0.074 - 0.324 


2.21 


2.42(03) 


0(6) 


2.4(5.3) 


3.4(2.1) 


0.074 - 0.324 


4.06 



Table 7: The parameters of the fits according to Eqs. (36) and (37). 
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Figure 8: The chiral condensate (ipip) at m q a = from Table 6 (filled circles) as a 
function of (3 C — j3. The dashed lines are fits to Eq. (36) using the points in the range 
0.05 < p e - (5 < 0.2 with (3 m from 0(2) (long dashes) and 0(6) (short dashes). The 
solid line is a fit to Eq. (37) with f3 m fixed to the 0(2)- value. 

with an ansatz including a correction-to-scaling term 

{W>)iP < (3c,m q a = 0) = B((3 C - + b x {[3 c - (3) e ] , (37) 

and all available points apart from the first one. In Fig. 8 we show the result for 
f3 m fixed to the 0(2)-value as a solid line. Here, the correction exponent is rather 
large, 9 ~ 2. A similar fit with /3 m (0(6)) is again worse and leads to large errors in 
the parameters. A fit to Eq. (37), where also (3 m is free, is very unstable and leads 
to a negative (3 m . Details of our fits are listed in Table 7. 
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4.3 Scaling behaviour of the chiral condensate 

In the last subsections we have seen that the chiral condensate data in the close 
vicinity of the chiral phase transition are compatible with 0(2) and, to a lesser 
degree, also 0(6) critical behaviour. In the following we want to compare the uni- 
versal scaling function for aQCD, / = (ip^m^ 1 ^ 5 , to those of the spin models, 
fa = Mhr l l & . In Fig. 9 we show the data for (ipip)m' 1 ^ 5 as a function of the scaling 
variable z = fir/m 1 ^ with the critical exponents for the 0(2) (a) and 0(6) (b) mod- 
els from Table 1. The used normalizations (m q a) and (3q are shown in Table 8. The 
mass normalization was calculated from the amplitude d c determined in subsection 
4.1. The normalization fl can be deduced from the amplitude B computed in the 
last subsection. Alternatively, one may fix the derivative of the scaling function at 
the critical point to the one of the corresponding 0(N) model. This is what we have 
done for the plots. The derivative is given by 

/'(O) = <£A)d c (1//3m ~ 1) • (38) 

The derivatives of the 0(2) and 0(6) scaling functions can be found from the respec- 
tive parametrizations in Refs. [23, 27]. Choosing this possibility we obtain B = 2.06 
for 0(2) and B = 2.27 for 0(6). That is, for 0(2) we find a 5-value which is in 
agreement with the 0(2) and free fit results from Table 7, whereas the 0(6) fits 
lead to somewhat higher 5-values. In Fig. 9 we have plotted only the data from 
the range 5.3 < (3 < 5.9, at f3 — 5.3 only the data for m q a < 0.03. All other data 
are either distorted by the deconfinement transition or too far outside the critical 
region. 

We observe in Fig. 9 (a) that the data scale very well with 0(2) exponents for 
small \z\ and agree there also with the 0(2) scaling function fa{ z )- For decreasing 
z < the data for different (3 start to deviate from each other. A similar effect was 
found for the 0(2) data [23] and is characteristic for corrections to scaling. At z > 
the data scale well apart from the data for m q a = 0.005. This may be due to the 
presence of finite size effects, though the few data with low statistics which we have 
from the larger lattices coincide inside the error bars with the 8 3 x 4 data. In Fig. 
9 (b), where we show the data using 0(6) parameters, the situation is somewhat 
different. Here, the data scale very well apart from the neighbourhood of the critical 





(m q a) 


A) 


d c 


dl 


&(0) 


B 


0(2) 


0.0093 


0.126 


2.66 


-0.366 


-0.286 


2.06 


0(6) 


0.0091 


0.145 


2.67 


-0.605 


-0.332 


2.27 



Table 8: The normalization constants and the corresponding amplitudes used in the 
scaling plots. 



20 




Figure 9: The scaling function f(z) = (4>?p) /m 1 / 5 versus the scaling variable z = 
(3 r / m 1//A with data from the 8 3 x 4 lattice for 0(2) (a) and 0(6) (b) critical exponents 
and normalizations from Table 8. The solid curves are the universal scaling functions 
from the respective 0(N) spin models. The data points are connected by straight 
lines to guide the eye, the numbers indicate the /3-values of the data. 
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Figure 10: Comparison of the scaling behaviour of {ipip) /m 1 ^ for 0(2) (left) and 
0(6) (right) exponents and normalizations in the vicinity of the critical point. The 
notation is the same as in the previous figure. For clarity the data for m q a = 0.005 
have been omitted. The black solid lines are the respective O(N) scaling functions. 



point where the data spread somewhat and deviate from the 0(6) scaling function, 
except very close to the normalization point z = 0. This is in contrast to the 0(2) 
case and can be seen in more detail in Fig. 10, which shows the scaling behaviour 
in the close vicinity of the critical point. 

As mentioned in the introduction there are first results [19] for the universality 
class which is expected to govern the chiral transition in the continuum limit. The 
universal scaling function of the class is still unknown, but estimates of two critical 
exponents, namely r\ ~ 0.2 , v ~ 1.1 are known. The errors are still large, e.g. 
r\ = 0.2(1) [33]. Due to the relation 

5 = ^1 (39) 

for d = 3 this implies 5 = 4.0(4) , and from the considerations in subsection 4.1 and 
the numbers in Table 3 this entails a shift in the location of the critical point. Though 
we have no scaling function to compare with we can still perform a test on the scaling 
of the data assuming critical exponents compatible with the estimates of Ref. [19]. 
We have varied the exponents correspondingly and calculated each time the adequate 
critical point. In Fig. 11 we show the best scaling result for ('0'0)/(m 3 a) 1//,5 as a 
function of ((3 - /3 c )/(m q a)^ A which is obtained for S = 4A,u= 1.0 (1/A = 0.409) 
and {3 C = 5.64. We see from Fig. 11 that the data still spread considerably around 
the critical point and do not indicate a definite scaling function. 
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Figure 11: Scaling test for {ipip) / (niga) 1 ^ with critical exponents from [19]. The 
data from the 8 3 x 4 lattice are plotted as a function of {(3 — /9 c )/(m q a) 1//A . The 
notation is the same as in Fig. 9. 

5 Summary and Conclusions 



In this paper we have studied QCD with two flavours of staggered Dirac fermions 
in the adjoint representation at finite temperature. Our main objective was the 
investigation of the chiral transition and the determination of the associated univer- 
sality class. For ordinary QCD with two staggered fermions this issue is still under 
discussion and may be difficult to decide because the deconfinement transition coin- 
cides with the chiral transition. In contrast to QCD the corresponding transitions of 
aQCD occur at two distinct critical temperatures and therefore allow to study the 
transitions separately. A major prerequisite of any scaling test on critical behaviour 
is the exact location of the critical point. In order to achieve that goal and also 
to obtain sufficient data for the tests we have extended the existing data base for 
aQCD supplied by the work of Karsch and Lutgemeier [16]. 

The critical point of the chiral transition was then determined using several 
strategies. Whereas the peak positions of the chiral susceptibility provide only first 
hints on a f3 c - value below 5.79 - the result of Ref. [16] - the evaluation of the mass 
dependence of (ipip) leads definitely to an interval 5.6 < (3 C < 5.7 and a value 
0.19 < 1/8 < 0.28 for the citical exponent of (ipip) on the critical line. The known 
values of 5 for the O(N) spin models are all close 4.8, that is 1/5 ~ 0.21. The 
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Figure 12: Three-dimensional plot of the chiral condensate (ipip) for all /3-values und 
m q a < 0.04. The data are from the 8 3 x 4 lattice. The dark blue lines connect the 
same masses, the green ones show fits to the critical behaviour on the coexistence 
line and the critical line. The red lines are fits of the mass dependencies at fixed 
/3 and the light blue lines show the expected behaviour in regions where we do not 
have enough data. The dotted lines show the hights and indicate the /9-values. 

value 5 = 4.0(4) of the new universality class proposed by Ref. [19] is however also 
compatible with the found 1/5-interval. In the following we improved on our (3 C - 
estimate by including the next term in the scaling ansatz, using the exponents of 
the 0(2) and 0(6) model and later also the ones from Ref. [19]. As it turned out 
both spin models lead to the same rather accurate value (3 C = 5.624(2). This value 
and the respective value (3d = 5.236(3) [16] of the deconfinement transition yield, 
via the two- loop beta function, the ratio T c /T d = 7.8(2). 

A further aim of our paper was the verification of the Goldstone effect which is 
expected to control the mass dependence of the chiral condensate below the chiral 
transition temperature. In the three-dimensional O(N) models the effect is even 
detectable in the critical region close to T c due to its influence on the universal 
scaling function fc(z). We have evaluated the chiral condensate data at fixed (3 < j3 c 
in the same manner as it was done for the magnetization of the three-dimensional 
0(N) models and find a completely analogous behaviour, that is aQCD acts as an 
effective three-dimensional O(N) spin model in its chirally broken phase. However, 
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close to the deconfinement transition this behaviour is disturbed and (ipip) develops 
a gap. Below (3 d) at (3 — 5.1, (ipip) becomes a smooth, though only slightly varying 
function of m q a, which could as well be interpreted as the behaviour of an effective 
four-dimensional theory. 

Using the obtained values for the critical point and the critical amplitudes for 
the normalization of the scales we have carried out scaling tests and comparisons to 
the scaling functions of the 0(2) and 0(6) models. With the 0(2) critical exponents 
the data are scaling convincingly in the direct neighbourhood of the critical point, 
that is for small values of \z\ and they coincide there with the 0(2) scaling function. 
For larger values of (—z) the data show - like in the original non-linear a-model - 
increasing deviations, which are typical for corrections to scaling. Using the 0(6) 
parameters for the scaling test seems at first sight to lead to scaling of the data in a 
larger region of z. A closer inspection in the vicinity of the critical point, where the 
0(2) test works, shows however, that the data spread there and deviate somewhat 
from the 0(6) scaling function. This observation is confirmed by the calculation of 
the amplitude B from the derivative of the scaling function: it is in better agreement 
with the value found from the Goldstone effect for 0(2) than for 0(6). We have 
also carried out scaling tests with the critical exponents calculated in Ref. [19] for 
the universality class of the continuum limit. Our lattice data do not yet support 
the proposed class. 

In conclusion, our findings have led to a consistent picture of the behaviour of 
aQCD from below the deconfinement transition to well above the chiral transition. 
We have found evidence for a second order chiral transition at (3 C = 5.624(2), com- 
patible with the 3d 0(2) universality class. The expected 3d Goldstone effect in 
the intermediate phase (3 d < (3 < (3 C , which was already noted in Ref. [16], was 
fully confirmed. The same applies to the first order deconfinement transition at 
(3 d = 5.236(3). Our results are visualized in Fig. 12, where we show the chiral con- 
densate data at small mass values together with fits based on our understanding of 
the behaviour of aQCD at finite temperature. 
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